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1 Introduction 

The study of molecular vibrational spectra [jj] requires theoretical models in order to analyze and interpret 
the measurements. These models range from simple parametrizations of the energy levels, such as the 
Dunham expansion to ab initio calculations, where solutions of the Schrodinger equation in different 
approximations are sought || ^, ^| . In general, the latter involve the use of internal coordinates and the 
evaluation of force field constants associated to derivatives at the potential minima. While this method 
can be reliably applied to small molecules 0, it quickly becomes a formidable problem in the case of 
larger molecules, due to the size of their configuration spaces. New calculational tools to describe complex 
molecules are thus needed. 

In 1981 an algebraic approach was proposed to describe the roto-vibrational structure of diatomic 
molecules |§| , subsequently extended to linear tri- and four- atomic molecules Q and certain non- linear 
triatomic molecules fic|| . Although these were encouraging results, the model could not be extended 
to polyatomic molecules, due to the impossibility of incorporating the underlying discrete symmetries. 
This difficulty could be surmounted by treating the vibrational degrees of freedom separately from the 
rotations. In 1984 Van Roosmalen et al. proposed a U(2) based model to describe the stretching 
vibrational modes in ABA molecules |Tl|] which was later extended to describe the stretching vibrations 
of polyatomic molecules such as octahedral and benzene like molecules |l2| . Recently the bending modes 
have also been included in the framework, which was subsequently applied to describe C2„-triatomic 
molecules |L3| and the lower excitations of tetrahedral molecules |l4|], using a scheme which combines 
Lie-algebraic and point group methods. In a different approach, it has also been suggested |D| to use a 
U(k + 1) model for the k = 3n — 3 rotational and vibrational degrees of freedom of a n-atomic molecule. 
This model has the advantage that it incorporates all rotations and vibrations and takes into account the 
relevant point group symmetry, but for larger molecules the number of possible interactions and the size 
of the Hamiltonian matrices increase very rapidly, making it impractical to apply. 

Although the algebraic formulations have proved useful, several problems remained, most important 
of which is the absence of a clear connection to traditional methods. On the other hand, a related 
problem is the lack of a systematic procedure to construct all physically meaningful interactions in the 
algebraic space. In this paper we show that both these issues can be resolved by considering a symmetry- 
adapted version of the U(2) algebraic model for the analysis of molecular vibrational spectra. In this 
approach it is possible to construct algebraic operators with well defined physical meaning, in particular 
interactions fundamental for the description of the degenerate modes present in systems exhibiting high 
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degree of symmetry. The procedure to construct them takes full advantage of the discrete symmetry of 
the molecule and gives rise to all possible terms in a systematic fashion, providing a clear-cut connection 
between the algebraic scheme and the traditional analyses based on internal coordinates, which correspond 
to the harmonic limit of the model |l6| . 

As a test of the symmetry-adapted approach we discuss an application to three 2?3^-triatomic molecu- 
lar systems, namely H^, Be3 and Naj" , and to two tetrahedral molecules, the Be4 cluster and the methane 
molecule. Since small molecules can in general be well described by means of ab initio calculations Jl7| , |T?| , 
we emphasize the basic purpose of this work. We establish an exact correspondence between configura- 
tion space and algebraic interactions by studying the harmonic limit of the U(2) algebra. This general 
procedure not only allows to derive algebraic interactions from interactions in configuration space, but 
can also be applied to cases for which no configuration space interactions are available. The 2? 3 ^-triatomic 
molecules constitute the simplest systems where degenerate modes appear and where the new interactions 
in the model become significant. In the case of Be4 we present a comparison with ab initio calculations, 
while for CH4 we present a detailed comparison with experiment. The application of these techniques to 
other molecules, as well as a more complete presentation can be found in references jl6| |l9[ f2l| . 



2 The U(2) vibron model 

The model is based on the isomorphism of the U (2) Lie algebra and the one dimensional Morse oscillator 

+ D{e- 2x ' d - 2e~ x/d ) , (1) 
. In order to see how this isomorphism 
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whose eigenstates £ can be associated with U{2) D SO{2) states 
comes about, consider the radial equation 
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r dr dr 



(r) = (N + l)0(r) 



which corresponds to a two-dimensional harmonic oscillator (in units where h = fi = e 

By carrying out a change of variable 
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1) associated 



to a U{2) symmetry algebra 
Eq. (2) transforms into 
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This can be identified with Eq. (1) after defining x = pd and multiplying by h 2 /2/id 2 , provided that 
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where we have defined m = a/2. In the framework of the U{2) algebra, the operator N corresponds to the 
total number of bosons and is fixed by the potential shape according to Eq. (4) , while m, the eigenvalue of 
the 5*0(2) generator J 2 , takes the values m = ±N/2, ±(N— 2)/2, .... The Morse spectrum is reproduced 
twice and consequently for these applications the m- values must be restricted to be positive. In terms of 
the £7(2) algebra, it is clear from Eqs. (3-5) that the Morse Hamiltonian has the algebraic realization 
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In addition, the U(2) algebra includes the raising and lowering operators J+ and </_, which connect 
different energy states, while the angular momentum operator is given by J 2 = N(N + 2)/4, as can be 
readily shown. 

The Morse Hamiltonian of Eq. (6) can be rewritten in the more convenient form 

H' = H + A^- = ^[{J + J_ + J_J + )-N] , (7) 

where we have used the relation J 2 = J 2 — ( J + J_ + J_ J+)/2 and added a constant term AN 2 / '4 in order 
to place the ground state at zero energy. The parameters N and A are related to the usual harmonic 
and anharmonic constants u> e and x e uj e used in spectroscopy. To obtain this relation it is convenient to 
introduce the quantum number 

N 

v = — -m, (8) 

which corresponds to the number of quanta in the oscillator. In terms of v, the corresponding energy 
expression takes the form 

N 2 A 

E' = -A(m 2 -—) = — 2 {N + 1/2) + A(N + l)(v + 1/2) - A(v + 1/2) 2 , (9) 
from which we immediately obtain 

cu e = A(N + 1) , 
x e u) e = A . (10) 

Thus, in a diatomic molecule the parameters A and iV can be determined by the spectroscopic constants 
Lo e and x e uj e . 

We now consider the U%{2) D SUi(2) D SOi{2) algebra, which is generated by the set {G,} = 
{Ni, J+s, J-.ii Jo,i}, satisfying the commutation relations 

[Ja,i,J±,i] — ± J±,i j [J+,i,J-,i] = 2Jo,j , [Ni, J M) j] = 0, (11) 

with fj, = ±,0. As mentioned before, for the symmetric irreducible representation [Ni, 0] of Ui(2) one can 
show that the Casimir operator is given by J 2 = Ni(Ni + 2)/4§, from which follows the identification 
ji = N/2. The SOi(2) label is denoted by rm. 



3 The Be4 cluster 

As a specific example, we consider the Be4 cluster, which has a tetrahedral shape. V^h molecules can 
be similarly treated. In the Be4 case there are six Ui{2) algebras involved (i = 1, . . . , 6). In the present 
approach each relevant interatomic interaction is associated with a Ui(2) algebra. The operators in the 
model are expressed in terms of the generators of these algebras, and the symmetry requirements of the 
tetrahedral group Td can be readily imposed |TJ, The local operators {Gi} acting on bond i can be 
projected to any of the fundamental irreps T = A\, E and Ft- Using the generators we obtain the 
T d tensors 

6 

where /i = ±, and 7 denotes the component of Y. The explicit expressions are given by 
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^(4,3-4,5) ■ (13) 



The Hamiltonian operator can be constructed by repeated couplings of these tensors to a total symmetry 
A\ , since it must commute with all operations in 7^ pif ■ 

All calculations are carried out in a symmetry-adapted basis, which is projected from the local basis 

t/i(2) <g> ••• g U 6 {2) D 5d(2) $ ••• ® 506(2) D 50(2) 

I II II (14) 

I [iV a ] , ... , [N 6 ] ; «i , . . . , « 6 ; K ) 

in which each anharmonic oscillator is well defined. By symmetry considerations, Ni = N for the six 
oscillators, Uj = Ni/2 — rrij denotes the number of quanta in bond i and V = 52, is the total number of 
quanta. The local basis states for each oscillator are usually written as \Ni,Vi), where V{ = (Ni — 2mi)/2 = 
0, 1, . . . [Ni/2] denotes the number of oscillator quanta in the t-th oscillator. The states with one quantum 
V = 1 are denoted by | i ) with Vi = 1 and Vj^i — 0. Using the same projection technique as for the 
generators ([13]), we find the six fundamental modes 

iX> = XX<i*>- (15) 
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The expansion coefficients are the same as in Eq. (|l3j ) . The states with a higher number of quanta 
can be constructed using the Clebsch-Gordan coefficients of 7^ |24| . Since all operators are expressed 
in terms of powers of the Ui(2) generators, their matrix elements can be easily evaluated in closed form. 
The symmetry-adapted operators of Eq. ( |l3| ) and symmetry-adapted basis states are the building blocks 
of the model. 

We now proceed to expicitly construct the Be4 Hamiltonian. For interactions that are at most 
quadratic in the generators the procedure yields 

H = uji H Al + UJ2 He + ^3 Wf 2 + a 2 V E + a 3 Vf 2 , (16) 

with 
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Vr = ~£*W„. (17) 



Note that we have not included Va x in H , since the combination 
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£(Wr + Vr) = -Lj2 Ni (N i+ 2), (18) 



r 1=1 



is a constant 3(iV + 2)/2 . The five interaction terms in Eq. ( ITq ) correspond to linear combinations of 
the Casimir operators of However, for a good description of the vibrational energies of Be4 it is 
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necessary to include interactions which are related to the vibrational angular momenta associated with 
the degenerate modes E and F2. These kind of terms is absent in the former versions of the model 
]l2| , |l4[ ]. We now proceed to show how they can be obtained in the present formalism. In configuration 
space the vibrational angular momentum operator for the E mode is given by |25| 

< 19 > 

where q^ and qf are the normal coordinates associated to the E mode. This relation can be transformed 
to the algebraic space by means of the harmonic oscillator operators 

v 6 -7i(^)' (20) 



to obtain 
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(ftftftf-ftftftf) • (21) 



Here = J^i a y % &») with a similar form for b^', while the can be read from Eqs. (12,13). In order 
to find the algebraic expression for l A2 we first introduce a scale transformation 

b\ ee J_,,;/0Vi , h ee S+^/y/N- . (22) 

The relevant commutator can then be expressed as 

Ml] = ^-[J + ,i,J-,i] = ^2J , 4 = , (23) 

where 

Vi = — - J ,i ■ (24) 

The other two commutation relations of Eq. (11) are not modified by the scale transformation of Eq. (22). 
In the harmonic limit, which is defined by Ni — ► 00, Eq. (23) reduces to the standard boson commutator 
= 1. This limit corresponds to a contraction of SU(2) to the Weyl algebra and can be used to 
obtain a geometric interpretation of algebraic operators in terms of those in configuration space. In the 
opposite sense, Eq. ( p2] ) provides a procedure to construct the anharmonic representation of harmonic 
operators through the correspondence fej — ► b\ = J-.i/VNi an d &i — * b% = J+^/\/Wi. Applying this 
method to the vibrational angular momentum we find 

i M = -jj(t-A f +*- f -,2 f *i) ■ ( 25 ) 

For the vibrational angular momentum l^ 1 associated with the F2 mode we find a similar expression. 
The corresponding interactions are 

Hi = g22i A2 i A2 + 9 33 J2 I i 1 ^ 1 ■ ( 26 ) 

7 

With this method we obtain an algebraic realization of arbitrary configuration space interactions. As a 
simple example, a one-dimensional harmonic oscillator Hamiltonian Hi = (b\bi + bib\)/2, transforms into 

^(J_,, J + ,i + J+,iJ-,i) = ^{Ji ~ Jli) = Vi + 1/2- ^ , (27) 
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where in the last step we used Eq. (24). The spectrum of Eq. (27) has an anharmonic correction, analogous 
to the quadratic term in the Morse potential spectrum. We are thus substituting harmonic oscillators by 
Morse oscillators. 

A more interesting application is to use our model to fit the spectroscopic data of several polyatomic 
molecules. In the case of Be4 the energy spectrum was analyzed by ab initio methods in [17]], where 
force-field constants corresponding to an expansion of the potential up to fourth order in the normal 
coordinates and momenta were evaluated. We have generated the ab initio spectrum up to three quanta 
using the analysis in pq ]. For the algebraic Hamiltonian we take |L6| 

H = uji H Al + oj 2 7~Le + <^3 Wj? 2 + X12 (ji-At T~Le^ + X 13 (Ha^f^^ + -^33 (t^f 2 ) 

+533 Z T ^ + *33 033 + t 2 3 623 , (28) 

7 

The terms O33 and O23 represent the algebraic form of the corresponding interactions in [ p5| which are 
responsible for the splitting of the vibrational levels in the {v\,v™ ,v\) — (0,0°, 2 2 ) and the (0, l 1 , l 1 ) 
overtones Jig ]. 

In Table I we show the the results of a least-square fit to the vibrational energies of Be4 with the 
Hamiltonian of Eq. (28). The r.m.s. deviation obtained is 2.6 cm -1 , which can be considered of spectro- 
scopic quality. We point out that in (25[ ^6) several higher order interactions are present which we have 
neglected. Since our model can be put into a one to one correspondence with the configuration space 
calculations, it is in fact possible to improve the accuracy of the fit considerably, but we have used a 
simpler Hamiltonian than the one of |2^, ^6). When no ab initio calculations are available (or feasible) 
the present approach can be used empirically, achieving increasingly good fits by the inclusion of higher 
order interactions 

We note that the Be4 Hamiltonian of Eq. (28) preserves the total number of quanta V. This is a good 
approximation for this case according to the analysis of [^5[ |2f|, but it is known that Fermi resonances 
can occur for certain molecules when the fundamental mode frequencies are such that (V, V') states with 
V 7^ V are close in energy. These interactions can be introduced in the Hamiltonian but the size of the 
energy matrices grows very rapidly, so the best way to deal with this problem is through perturbation 
theory. 



4 Uzh triatomic molecules 

For 2? 3 h molecules we follow a similar procedure, namely, we construct the symmetry-adapted op- 
erators and states analogous to Eq. (13,15) and carry out the building up procedure to construct the 
Hamiltonian and states with a higher number of quanta with the appropriate projection operators and 
Clebsch-Gordan coefficients pjjfl . 

In Table II we present the fits to the spectra of Be3, Najj" and H^~ up to three quanta. While remarkably 
accurate descriptions of the first two molecules can be achieved using a four-parameter Hamiltonian, we 
were forced to include four additional higher order terms in the H^~ Hamiltonian in order to properly 
describe this molecule. This is in accordance with the work of Carter and Meyer who were forced 
to include twice as many terms in the potential energy surface for H^~ than for the NaJ molecule. The 
H^ ion is a very "soft" molecule which, due to the light mass of its atomic constituents carries out large 
amplitude oscillations from its equilibrium positions |Tq] . 



5 The methane molecule 

We now turn our attention to the CH4 molecule, for which we shall make a detailed description. For 
methane we have four U (2) algebras corresponding to the C-H interactions and six more representing the 
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H-H couplings. The assignments and the choice of the Cartesian coordinate system are the same as in 
The molecular dynamical group is then given by the product 

G = E/i(2)® f/ 2 (2)®...<g) f/i (2) . (29) 

The labeling is such that i = 1, ...,4 correspond to the C-H couplings while the other values of i 
are associated with H-H interactions. Consequently there are two different boson numbers, N s for the 
C-H couplings and Nb for the H-H couplings, which correspond to the stretching and bending modes, 
respectively. The tctrahedral symmetry of methane is taken into account by projecting the local operators 
{Gj}, which act on bond i, on the irreducible representations T of the tetrahedral group 7^- The explicit 
expressions for the 7^ tensors for the stretching modes are 
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(30) 



while for the bending modes we have 
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T.,r) — [ e/u.fi J It. 6 Ji 



T F2 -» ~ — f T f 



M>3 



= ^=(4,7-4,9) • (31) 

As before, the algebraic Hamiltonian can be constructed by repeated couplings of these tensors to a total 
symmetry A\. 

The methane molecule has nine vibrational degrees of freedom. Four of them correspond to the 
fundamental stretching modes {A\ © Fg) and the other five to the fundamental bending modes (E © F2) 
^?J. The projected tensors of Eqs. ( |30| ) and ( fjlf ) correspond to ten degrees of freedom, four of which 
(A\ © F2) are related to stretching modes and six (Ai © E © F2) to the bendings. Consequently we 
can identify the tensor T^' 6 as the operator associated to a spurious mode. This identification makes 

it possible to eliminate the spurious states exactly. This is achieved by (i) ignoring the T^ b tensor 
in the construction of the Hamiltonian, and (ii) diagonalizing this Hamiltonian in a symmetry-adapted 
basis from which the spurious mode has been removed following the procedure of flij . We note that 
the condition on the Hamiltonian that was used in to exclude the spurious contributions, does not 
automatically hold for states with higher number of quanta. 

According to the above procedure, we now construct the 7^ invariant interactions that are at most 
quadratic in the generators and conserve the total number of quanta 
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^ = atE^^. (32) 
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Here r = A±, F2 for the stretching vibrations x = s and r = E, F2 for the bending vibrations x — b. In 
addition there are two stretching-bending interactions 

»- - < 33 > 

The zeroth order vibrational Hamiltonian is now written as 

H a = Ldi Ha 1>3 + HE b + W3 Hf 2 , s + w 4 HiT 2 + UJ34 Hsb 

+a 2 V Eb + a 3 Vf 2 , s + a 4 VF 2 . b + a 34 V sb . (34) 
The interaction Va x , has not been included since, in analogy to Eq. (18), the combination 

, 4 

r s i=i 

corresponds to a constant N s + 2. A similar relation holds for the bending interactions, but in this case 
the interaction VUi.i, has already been excluded in order to remove the spurious A\ bending mode. The 
subscripts of the parameters correspond to the (2/1, v l 2 2 , v 1 ^ , v 1 ^) labeling of a set of basis states for the 
vibrational levels of CH4. Here v\ 1 1/2, ^3 and ^4 denote the number of quanta in the Ai_ s , Ef,, i^.s and 
F2.b modes, respectively. The labels U are related to the vibrational angular momentum associated with 
degenerate vibrations. The allowed values are k = i/j, z/j — 2, . . . , 1 or for vi odd or even | p7| . 

In the harmonic limit the interactions of Eqs. (32) and (33) again attain a particularly simple form, 
which can be directly related to configuration space interactions. This limit is obtained, as before, by 
rescaling J +i i and by y/Wl and taking A/, — > 00, so that 

lim —= = bi , 



lim — == = b 



lim —[j +l J_ l ] = lim ^1 = 1. (3 6 ) 

where the operators bi and W satisfy the standard boson commutation relation [6^, &t] = Sij. Applying 
the harmonic limit to the interactions of Eqs. (32) and (33) we obtain 

7 

lim Vr x = , 

N x —>oo 
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lim V sb = . (37) 



Here the operators t are given in terms of the local boson operators b\ through the coefficients a 1 



given in Eqs. (30, 3f ) 



s,7 'bj-y 



6T-t = ga^&t, (38) 
i=l 

with a similar relation for the annihilation operators. From Eq. (37) the physical interpretation of 
the interactions is immediate. The 7i.r x terms represent the anharmonic counterpart of the harmonic 
interactions, while the Vr^ terms are purely anharmonic contributions which vanish in the harmonic 
limit. In an application to the ozone molecule it was found that these terms can account for the strong 
anharmonicities and incorporate the effect of Darling-Dennison type couplings p0| . 

The zeroth order Hamiltonian of Eq. (34) is not sufficient to obtain a high-quality fit of the vibrations 
of methane. Several physically meaningful interaction terms that are essential for such a fit are not 
present in Eq. (34). They arise in our model as higher order interactions. It is an advantage of the 
symmetry-adapted model that the various interaction terms have a direct physical interpretation and a 
specific action on the various modes [|l6) . Hence the addition of higher order terms and anharmonicities 
can be done in a systematic way. For the study of the vibrational excitations of methane we use the 7~d 
invariant Hamiltonian |2l[| 

H = ui Ha 1iS + U2 "He* + k>3 Tifa + w 4 Hf 2 , 6 + a 3 V F2 s 

+Xn (Ha 1iS j + X 22 (li.E b ^J + X 33 \Hf 2 ,s^ +X ii {ii F2b 

+X12 {Ha 1-s 7~LE b ^J + X u (Ha 1:S ^F 2 , b ) 

+X 23 {^E b 7^F 2 , S ) + X 2 4 (ji-E,, T~Lf 2 ^ + X 3i (Hf2 }3 T^F; 

+.922 (l A >) 2 + <? 33 £ I* I* + 544 E %n + ^ E 1 

777 

+t 33 o ss + t iA 6 bb + t 3i 6 sb + 1 23 6 2s + t 24 6 2b . (39) 

The interpretation of the and a 3 terms follows from Eq. (37). The Xij terms are quadratic in the 
operators H-r^ and hence represent anharmonic vibrational interactions. The gtj terms are related to 
the vibrational angular momenta associated with the degenerate vibrations. As mentioned before, these 
interactions, which are fundamental to describe molecular systems with a high degree of symmetry, are 
absent in previous versions of the vibron model in which the interaction terms are expressed in terms of 
Casimir operators and products thereof jl^, |l4| . They give rise to a splitting of vibrational levels with 
the same values of (yi, v 2 , v 3 , u 4 ) but with different l 2 , l 3 and/or I4. Their algebraic realization is given 
by 

l M = -i V2—[f Eb x f+ b ) A2 , 
N b 

1* = +l V2^-[f F ^xf^ . (40) 
The square brackets in Eq. (40) denote the tensor coupling under the point group 7^ 

l^'xf 2 ]; = E ^'^^571,72,7)^^, (41) 

7l>72 

where the expansion coefficients are the Clebsch-Gordan coefficients for 7^ ^4j. In the harmonic limit 
the expectation value of the diagonal terms in Eq. (39) leads to the familiar Dunham expansion [^7) 

E~< :: ''< ' • XX x n ^ + + y) + EE an Uh ■ (42) 

i j>i i j>i i 
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Here di is the degeneracy of the vibration. The t%j terms in Eq. (39) give rise to further splittings of the 
vibrational levels (v\, v 2 , ^3, ^4) into its possible sublevels. They can be expressed in terms of the tensor 
operators of Eqs. (30) and (31) as 

6 xy = ( 6]T[T x f F2 '»]Z [Tf ■» x - 4^[f x f F2 X 2 \T+- y x t?~]** J , 

3 a \ 7 7 / 

a 2;c = ^^fa^T^xTf 3 '-]*^ . (43) 

1 \ 7 7 / 

In the harmonic limit the t%j terms have the same interpretation as in p5| . The O ss , Obb and O s b terms 
give rise to a splitting of the E and F 2 vibrations belonging to the {v\ , v£ , v£ , v£ ) = (0, 0°, 2 2 , 0°), 
(0,0°,0°,2 2 ) and (0, 0°, l 1 , l 1 ) levels, respectively. Similarly the 6 2s and 6 2b terms split the Ft and F 2 
vibrations belonging to the (0, l 1 , l 1 , 0°) and (0, l 1 , 0°, l 1 ) overtones, respectively. 

The Hamiltonian of Eq. (^) involves 23 interaction strengths and the two boson numbers, N s and 
Nb. The vibron number associated with the stretching vibrations is determined from the spectroscopic 
constants cu e and x e uj e for the CH molecule to be N s — 43 ^3|. The vibron number for the bending 
vibrations, which are far more harmonic than the stretching vibrations, is taken to be iV& = 150. We 
have carried out a least-square fit to the vibrational spectrum of methane including 44 experimental 
energies from p8 |. The values of the fitted parameters are presented in the second column of Table III 



(Fit 1). In Table IV we compare the results of our calculation with the experimentally observed energies. 
All predicted levels up to V = 3 quanta are included. We find an overall fit to the observed levels with 
a r.m.s. deviation which is an order of magnitude better than in previous studies. Whereas the r.m.s. 
deviations of standard vibron model studies in and |^9| are 12.16 and 11.61 cm -1 for 19 energies, we 
find a value of 1.16 cm -1 for 44 energies. 

In order to address the importance of the 0:3 term, which is completely anharmonic in origin and 
vanishes in the harmonic limit (see Eq. (J37|)), we have carried out a calculation without this term (Fit 2 



of Tabic [II). With one less interaction term the r.m.s. deviation increases from 1.16 to 4.49 cm -1 . This 
shows the importance of the 013 term to obtain an accurate description of the anharmonicities that are 
present in the data. Similarly, the third calculation (Fit 3) shows that, in a fit without the t%j terms, the 
r.m.s. deviation increases from 1.16 to 7.81 cm -1 . 

In addition, we have carried out a fit in the harmonic limit (N s , Nb — > 00). In this limit the a 3 
term vanishes and the algebraic Hamiltonian of Eq. |3^) reduces to the vibrational Hamiltonian of 0] , 
the harmonic frequencies u>i and anharmonic constants Xy, and tij having the same meaning. The 
r.m.s. deviation increases to 20.42 cm -1 (Fit 4). When we use the vibrational Hamiltonian of which 



contains one additional interaction (called S34 in Table III) the quality of the fit does not improve. In 
fact, the r.m.s. deviation increases to 20.90 cm -1 (Fit 5). The importance of the nondiagonal elements is 
demonstrated in a calculation (Fit 6), in which the data are fitted by the (diagonal) Dunham expansion 
of Eq. ([42|). For this case the r.m.s. deviation is 21.00 cm -1 , almost the same value as for the other 
two calculations in the harmonic limit. The small differences in the parameter values and the r.m.s. 
deviation of Fits 4 and 6 show that, unlike for finite N s and Nb (Fits 1 and 3), in the harmonic limit the 
nondiagonal contributions from the tij terms are not very important. 



A comparison between the parameter values and the r.m.s. deviations of Fits 1-6 in Table III shows 
that the 0:3 term and the anharmonic effects in the interaction terms of Eq. ( |39| ) can only be compensated 
for in part by the anharmonicity constants Xij. The r.m.s. deviation increases from 1.16 to 4.49 and 
20.42 cm -1 for Fits 1, 2 and 4, respectively. 



6 Summary and conclusions 

In summary, in this paper we have studied the vibrational excitations of several molecules in a symmetry- 
adapted algebraic model. In particular, for the methane molecule we find an overall fit to the 44 observed 
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levels with a r.m.s. deviation of 1.16 cm -1 , which can be considered of spectroscopic quality. We pointed 
out that for this calculation the term in Eq. (39) in combination with the anharmonic effects in the 
other interaction terms plays a crucial role in obtaining a fit of this quality. Purely anharmonic terms 
of this sort arise naturally in the symmetry-adapted algebraic model, but vanish in the harmonic limit. 
Physically, these contributions arise from the anharmonic character of the interatomic interactions, and 
seem to play an important role when dealing with molecular anharmonicities, especially at higher number 
of quanta. This conclusion is supported by our other applications of the symmetry-adapted model to 
the Be4 cluster |l(| and the Hjj" , Be3 and NaJ molecules |lj| , as well as our study of two isotopes of the 
ozone molecule |20) . 

These studies suggest that the symmetry-adapted algebraic model provides a numerically efficient tool 
to study molecular vibrations with high precision. The main difference with other methods is the use 
of symmetry-adapted tensors in the construction of the Hamiltonian. In this approach, the interactions 
can be constructed in a systematic way, each term has a direct physical interpretation, and spurious 
modes can be eliminated exactly. It will be important to further explore the scope and applicability of 
the present approach. A more extensive study of methane including rotation-vibration couplings, states 
with a higher number of quanta and transition intensities is in progress. 
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Table I: Fit to ab initio E(J calculations for Be4. The values of the parameters are U)i = 636, U2 = 453, uj 3 = 532, 
X33 = 44.276, Xia = 4.546, X 13 = -2.539, g 33 = -15.031, t 33 = -1.679 and t 23 = -1.175. The total number of 
bosons is TV = 44. The parameters and energies are given in cm - . 



V 




r 


Ab initio 


Fit 


V (vt,^,v l 3 ) 


r 


Ab initio 


Fit 








N -> 00 


TV = 44 






N -> 00 


N = 44 


1 


(1,0°, 0°) 


A x 


638.6 


637.0 


3 (1,0°, 2°) 


At 


2106.8 


2105.6 




(o,i\o°) 


E 


453.6 


455.0 


(1,0°, 2 2 ) 


E 


2000.1 


1999.8 




(O.OM 1 ) 


F-2 


681.9 


678.2 




F 2 


2056.8 


2052.8 


2 


(2,0°,0°) 




1271.0 


1269.2 


(0,3\0°) 


E 


1341.3 


1343.7 




(l,l\0°) 


E 


1087.1 


1087.0 


(0,3 3 ,0°) 


At 


1355.5 


1352.5 




(1,0m 1 ) 


F-2 


1312.6 


1308.3 




A 2 


1355.5 


1354.4 




(0,2°,0°) 


A t 


898.3 


901.4 


(o^-M 1 ) 


F 2 


1565.5 


1565.7 




(0,2 2 ,0°) 


E 


905.4 


906.1 




F 2 


1584.4 


1583.1 




(Q, l 1 , l 1 ) 


Fi 


1126.7 


1125.1 


(0,2m 1 ) 


Fi 


1578.5 


1578.0 






F 2 


1135.5 


1134.1 


(0,1\2°. 2 ) 


E 


1821.4 


1821.6 




(0,0°, 2°) 


Ai 


1484.0 


1483.0 




E 


1929.5 


1929.0 




(0,0°, 2 2 ) 


E 


1377.3 


1373.9 


(0,1\2 2 ) 


A 2 


1813.3 


1813.1 






F 2 


1434.1 


1429.6 




At 


1830.8 


1831.7 


3 


(3,0°,0°) 


A t 


1897.0 


1896.7 




F 2 


1874.4 


1873.2 




(2,1 1 ,0°) 


E 


1714.3 


1714.3 




Fi 


1883.2 


1883.0 




(2,0°, l 1 ) 


F-2 


1937.0 


1933.7 


(0,0°, S 1 ^ 3 ) 


F 2 


2136.5 


2134.2 




(1,2°,0°) 


At 


1526.6 


1529.2 




F 2 


2327.3 


2326.9 




(1,2 2 ,0°) 


E 


1533.7 


1532.8 


(0,0°, 3 3 ) 


Ft 


2199.8 


2197.1 




(1 l 1 l 1 ) 


Ft 


1752.2 


1749.7 




At 


2256.5 


2254.4 






F 2 


1761.0 


1759.8 
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Table II: Least-square energy fit for the vibrational excitations of Hj, Be3 and NaJ. The energy differences 
AE = Eth — E exp are given in cm -1 . 











Be 3 


Na+ 


V 




r 


AE 


AE 


AE 


1 


(O,! 1 ) 


E 


-1.55 


0.51 


0.93 




(1,0°) 




0.42 


0.02 


1.95 


2 


(0,2°) 


M 


7.48 


-0.74 


0.37 




(0,2 2 ) 


E 


-5.69 


0.17 


0.84 




(I,! 1 ) 


E 


-0.61 


0.82 


1.68 




(2,0°) 


At 


-0.11 


-0.04 


1.26 


3 


(0.3 1 ) 


E 


-4.46 


-2.05 


-1.19 




(0,3 3 ) 


A x 


3.18 


-1.23 


-0.34 




(0,3 3 ) 


A 2 


2.44 


0.61 


-0.33 




(1,2°) 


Ai 


0.66 


1.90 


-0.01 




(1,2 2 ) 


E 


-5.00 


-1.36 


0.34 




(2, I 1 ) 


E 


4.07 


0.79 


-0.19 




(3,0°) 


Ax 


-1.23 


-1.66 


-2.06 






r.m.s. 


5.84 


1.35 


1.33 



Parameters 8 4 4 
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Table III: Parameters in cm 1 obtained in the fits to the vibrational energies of CH4. 



Parameter 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


Fit 5 


Fit 6 


N s 


43 


43 


43 


00 


00 


00 


N b 


150 


150 


150 


00 


00 


00 


Wi 


2977.60 


2966.17 


2970.03 


2967.40 


2966.81 


2969.75 


U)2 


1554.83 


1549.96 


1550.88 


1558.38 


1558.51 


1548.30 


U> 3 


3076.45 


3076.41 


3079.33 


3081.34 


3082.00 


3060.48 


LO4 


1332.22 


1329.52 


1337.69 


1337.51 


1337.54 


1338.64 


a 3 


582.87 


- 


480.33 


- 


- 


- 


Xu 


3.69 


10.23 


6.06 


-21.30 


-21.19 


-21.59 


X22 


1.30 


1.32 


1.37 


-1.17 


-1.17 


-0.28 


X33 


5.43 


6.97 


5.34 


-10.79 


11.12 


-7.89 


X44 


-3.47 


-3.64 


-4.41 


-6.26 


-6.27 


-7.26 


X\2 


-3.60 


-0.94 


-1.47 


-3.39 


-3.28 


-2.80 


X\3 














X\4 


-2.86 


-0.49 


-2.30 


-3.10 


-3.00 


-4.48 


X23 


11.14 


-8.75 


-10.68 


-7.97 


-8.10 


-3.80 


X24 


1.00 


0.91 


2.03 


-5.37 


-5.37 


-4.74 


X34 


-5.60 


-3.97 


-6.50 


-3.46 


-3.50 


-1.21 


922 


-0.46 


-0.46 


-0.41 


0.37 


0.37 


-0.62 


933 


0.19 


-1.23 


0.25 


-4.35 


-4.22 


-5.49 


#44 


4.07 


4.11 


3.79 


4.98 


4.98 


5.30 


334 


-0.65 


-0.72 


-0.56 


-0.74 


-0.87 


-0.67 


^33 


0.40 


0.16 




-1.25 


-1.27 




t 4 4 


1.00 


1.00 




0.56 


0.56 




^34 


0.21 


0.24 




0.24 


0.25 




^23 


-0.39 


-0.39 




-0.39 


-0.39 




^24 


0.13 


0.13 




0.91 


0.91 




S34 










-0.11 




r.m.s. 


1.16 


4.49 


7.81 


20.42 


20.90 


21.00 
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Table IV: Fit to vibrational excitations of CH4. The values of the parameters are given in the second column of 

E ca l 



Table|IIIl Here AE 
in cm 



E eX p- The experimental energies are taken from 



The wave numbers are given 



-1 



r 


(Vi . I/O. Vl, VA } 


E 1 


E 

exp 


AE 


r 


(Vl . 1/9 . Vi. . VA ) 

V -Li ii) 05 


E 1 


E 

exp 


AE 


4 1 


(1000) 


2916.32 


2916.48 


-0.16 




(01 1 1 ) 


5844 98 






E 


(0100) 


1533.46 


1533.33 


0.13 




(1200) 


5974.81 






F% 


(0001) 


1309.86 


1310.76 


-0.90 




(1011) 


7147.49 








(0010) 


3018.09 


3019.49 


-1.40 




(0021) 
(2100) 


7303.38 
7315.60 






M 


(0002) 


2587.77 


2587.04 


0.73 




(0120) 


7479.48 








(0200) 


3063.66 


3063.65 


0.01 




(0120) 


7557.17 








(0011) 


4323.81 


4322.72 


1.09 




(1020) 


8833.05 








(2000) 


5790.13 


5790 


0.13 


Fi 


(0003) 


3920.46 


3920.50 


-0.04 




(0020) 


5966.57 


5968.1 


-1.53 




(0102) 


4128.38 


4128.57 


-0.19 


E 


(0002) 


2624.14 


2624.62 


-0.48 




(0201) 


4364.39 


4363.31 


1.08 




(0200) 


3065.22 


3065.14 


0.08 




(0012) 


5620.08 








(0011) 


4323.09 


4322.15 


0.94 




(0012) 


5630.76 








(1100) 


4446.41 


4446.41 


0.00 




(1101) 


5755.58 








(0020) 


6045.03 


6043.8 


1.23 




(0111) 


5829.79 






Ei 


(0101) 


2845.35 


2846.08 


-0.73 




(0111) 


5848.94 








(0011) 


4323.15 


4322.58 


0.57 




(0210) 


6061.57 








(0110) 


4537.57 


4537.57 


0.00 




(10111 


7147.53 






E 2 


(0002) 


2612.93 


2614.26 


-1.33 




(0021) 


7303.29 








(0101) 


2830.61 


2830.32 


0.29 




(0021) 


7343.21 








(1001) 


4223.46 


4223.46 


0.00 




(1110) 


7361.79 








(0011) 


4321.02 


4319.21 


1.81 




(0120) 


7518.70 








(0110) 


4543.76 


4543.76 


0.00 




(0030) 


8947.65 


8947.95 


-0.30 




(1010) 


5845.53 






F 2 


(0003) 


3871.29 


3870.49 


0.80 




(0020) 


6003.65 


6004.65 


-1.00 




(0003) 
(0102) 


3931.36 
4143.09 


3930.92 
4142.86 


0.44 
0.23 


4, 


(0003) 


3909.20 


3909.18 


0.02 




(0201) 


4349.01 


4348.77 


0.24 




(0102) 


4131.92 


4132.99 


-1.07 




(0201) 


4378.38 


4379.10 


-0.72 




(0300) 


4595.26 


4595.55 


-0.29 




(1002) 


5523.80 








(1002) 


5498.66 








(0012) 


5594.92 


5597.14 


-2.22 




(0012) 


5617.16 








(0012) 


5620.68 








(0111) 


5836.11 








(0012) 


5632.36 








(1200) 


5973.26 








(1101) 


5740.86 








(1011) 


7147.56 








(out) 


5830.28 








(0021) 


7300.85 








(0111) 


5848.46 








(0120) 


7562.91 








(0210) 


6054.58 








(3000) 


8583.81 








(0210) 


6067.03 








(1020) 


8727.97 








(2001) 


7094.16 








(0030) 


8975.64 


8975.34 


0.30 




(1011) 


7145.84 






A 2 


(0102) 


4161.52 


4161.87 


-0.35 




(0021) 


7266.11 








(0300) 


4595.28 


4595.32 


-0.04 




(0021) 


7303.38 








(0111) 


5844.61 








(0021) 


7344.87 








(0120) 


7550.53 








(1110) 


7365.83 






E 


(0102) 


4105.22 


4105.15 


0.07 




(0120) 


7514.67 








(0102) 


4152.15 


4151.22 


0.93 




(2010) 


8594.90 








(0300) 


4592.13 


4592.03 


0.10 




(1020) 


8786.05 








(1002) 


5535.04 








(0030) 


8907.91 


8906.78 


1.13 




(0012) 


5620.36 








(0030) 


9045.36 


9045.92 


-0.56 




(0111) 


5836.45 
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